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On the ability of various circular inspiral templates that incorporate radiation 
reaction effects at the second post-Newtonian order to capture inspiral gravitational 
waves from compact binaries having tiny orbital eccentricities 

Manuel TessmeiQ and Achamveedu GopakumaiQ 
Theoretisch-Physikalisches Institut, Friedrich-Schiller-Universitat Jena, Max-Wien-Platz 1, 07743 Jena, Germany 

(Dated: December 2, 2008) 

We probe the ability of various types of post-Newtonian(PN)-accurate circular templates to cap- 
ture inspiral gravitational-wave (GW) signals from compact binaries having tiny orbital eccentrici- 
ties. The GW signals are constructed by adapting the phasing formalism, available in T. Damour, 
A. Gopakumar, and B. R. Iyer, [Phys. Rev. D 70, 064028 (2004)]. And, we employ the orbital en- 
ergy and the time-eccentricity to describe the orbital evolution of eccentric binaries that incorporate 
the 3PN-accurate conservative and the 2PN-accurate reactive dynamics. Using the fitting factor 
estimates, relevant for the initial LIGO, we show that circular templates, based on the adiabatic 
TaylorTl, complete adiabatic TaylorTl and TaylorT4 approximants, that require the 2PN-accurate 
energy flux are unable to capture our GW signals from compact binaries having tiny residual orbital 
eccentricities. However, the 2PN-order circular inspiral templates based on the recently introduced 
TaylorEt approximant are found to be both effectual and faithful in capturing GWs from inspiralling 
compact binaries having moderate eccentricities and we provide physical explanations for our obser- 
vations. We conclude that further investigations involving the actual interferometric data may be 
required to probe the ability of the widely employed traditional 2PN order circular templates to cap- 
ture inspiral GWs from astrophysical compact binaries, like the plausible compact binary progenitor 
candidate for the GRB 070201, that are expected to contain tiny residual orbital eccentricities. 

PACS numbers: 04.30.Db, 04.25.Nx 04.80.Nn, 97.60.Jd, 95.55Ym 



I. INTRODUCTION 



Inspiral GWs from non-spinning stellar mass compact 
binaries are being searched in the data generated by the 
operational first generation ground-based laser interfer- 
ometric GW detectors The GW data analysts em- 
ploy the technique of matched filtering that demands 
theoretically modeled inspiral templates to extract astro- 
physical GW signals from the noisy interferometric data. 
The theoretical inspiral templates require the temporally 
evolving GW polarizations h + (t) and h x (t), and are usu- 
ally constructed by employing the PN approximation to 
general relativity. In the case of inspiralling compact bi- 
naries, usually modeled to consist of point masses, the 
PN approximation provides, for example, the orbital dy- 
namics as corrections to the Newtonian equations of mo- 
tion in terms of {v/c) 2 ~ Gm/c 2 r, where v, m and r are 
the characteristic orbital velocity, the total mass, and the 
typical orbital separation, respectively. 

The LIGO Scientific Collaboration (LSC) 0] employs 
various types of PN-accurate inspiral templates, detailed 
in Ref. [J], for analyzing the LIGO data for inspiral GWs. 
These templates are constructed by keeping only the 
Newtonian contributions to the amplitudes of h + (t) and 
h x (t), while their phase evolutions are PN-accurate, re- 
sulting in the so-called restricted PN waveforms. While 
searching for inspiral GW signals from non-spinning com- 
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pact binaries, it turned out that the LSC had invoked 
various types of restricted 2PN order waveforms 
These inspiral template families provide slightly differ- 
ent prescriptions for the GW phase evolution that incor- 
porate gravitational radiation reaction (RR) effects at 
the 2PN order, i.e., (v/c) corrections beyond the New- 
tonian (quadrupolar) RR order. These templates rep- 
resent gravitational radiation field generated by compact 
binaries inspiralling under the action of 2PN-accurate RR 
effects along circular orbits. The approximation of quasi- 
circularity, namely inspiral along a sequence of circular 
orbits, is justified due to the fact that the gravitational 
RR forces are highly efficient in circularizing orbits of 
compact binaries [6(. 

Interestingly, the ability of the 2PN-accuratc circular 
templates, available in Ref. to detect GWs from com- 
pact binaries in inspiralling eccentric orbits is not ex- 
plored in the literature and we speculate that it may be 
due to Ref. 0. With the help of GW signals associated 
with compact binaries moving in Newtonian-accurate ec- 
centric orbits and inspiralling under quadrupolar RR or- 
der, Ref. Q demonstrated that the Newtonian order 
circular templates are quite efficient in extracting even 
their mildly eccentric GW signals. However, we recently 
pointed out that the arguments presented in Ref. Q 
should not be used to support the ability of the tradi- 
tional PN-accurate circular templates to capture inspiral 
GWs from compact binaries having tiny orbital eccentric- 
ities Q . This is because in Ref. Q , we demonstrated that 
the three types of circular inspiral templates based on the 
adiabatic, complete adiabatic and gauge-dependent com- 
pletely non-adiabatic approximants, detailed in Ref. Q 
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and relevant for the circular inspiral under the quadrupo- 
lar RR order, are inefficient in capturing GWs from ec- 
centric compact binaries having 2.5PN-accurate orbital 
evolution, detailed in Ref. fioj . The present paper ex- 
tends the analysis presented in Ref. Q by incorporating 
RR effects to 2PN (relative) order, while constructing 
both the GW signals and the circular search templates. 

In this paper, we explore the ability of a number of 
PN-accurate circular templates, detailed in Refs. [J, 0, 
[Til . [l2l ] , that employ the radiation reaction effects to the 
relative 2PN accuracy to capture GW signals from com- 
pact binaries in inspiralling eccentric orbits, constructed 
by adapting the phasing formalism, detailed in Ref. [IcI ]. 
We construct GW signals by perturbing compact bi- 
naries moving in 3PN-accurate processing eccentric or- 
bits with RR effects that are 2PN-accurate, beyond the 
quadrupolar (Newtonian) order. Further, we employ the 
orbital binding energy and the time-eccentricity et, ap- 
pearing in the PN-accurate generalized quasi-Keplcrian 
parametrization for eccentric orbits [l3|, to characterize 
the inspiralling and precessing eccentric binaries. We ex- 
plore the ability of various 2PN order circular templates 
to capture GW signals from compact binaries in inspi- 
ralling eccentric orbits by computing the fitting factors 
(FF), described in Ref. [14j . In the present study, the 
various 2PN order inspiral templates are provided by the 
adiabatic and complete adiabatic Taylor Tl, the adia- 
batic Taylor T4 and the TaylorEt approximants [see 
Eqs. (EH, tU, OH) and (IBJ) ]. 

The organization of the paper is as follows. In Sec. [Til 
we briefly describe the phasing of GWs from com- 
pact binaries in inspiralling eccentric orbits, available in 
Ref. ■ After that we provide certain PN-accurate for- 
mulae and relations that are necessary to construct our 
GW signals with the help of the lengthy PN-accurate ex- 
pressions, available in Ref. Section [TTJ also contains 
formulae required to construct various types of circular 
templates. How and why we perform a number of FF 
computations are explained in Sec. IIIII along with our 
results. Our conclusions are presented in Sec. HVl 



II. PHASING FORMULAE REQUIRED TO 
MODEL GW SIGNALS FROM ECCENTRIC 
BINARIES AND VARIOUS CIRCULAR 
TEMPLATES 

This section provides brief descriptions and necessary 
formulae required to model GWs associated with non- 
spinning compact binaries inspiralling along eccentric 
and circular orbits. For the sake of simplicity, we only 
consider temporally evolving restricted h x (t), having 
Newtonian accurate amplitude and PN-accurate phase 
evolution. In the case of compact binaries in inspiralling 
eccentric orbits, as noted earlier, we construct h x (t) 
by adapting the GW phasing formalism, available in 
Ref. [1 01 ] . This approach provides an efficient way of mod- 
cling the PN-accurate orbital dynamics of eccentric bina- 



ries incorporating the three inherent time scales, namely, 
those associated with the radial motion (orbital period), 
advance of periastron, and radiation reaction, crucial to 
construct temporally evolving h x (t). At present, the de- 
tailed PN-accurate expressions, describing the orbital dy- 
namics of non-spinning compact binaries in 3PN-accuratc 
eccentric orbits perturbed by the relative 2PN-accurate 
RR effects, computed with the help of Refs. [1(1 EH, are 
available in Ref. [l5| . In the next subsection, we summa- 
rize what were detailed in Refs. [HEHEH and obtain the 
restricted h x {t) associated with non-spinning compact 
binaries in precessing and inspiralling eccentric orbits in 
harmonic gauge. Due to the availability of the lengthy 
PN-accurate expressions in Ref. [l5| that are useful to 
construct our fiducial GW signals, wc only provide certain 
PN-accurate relations that are crucial to implement the 
the present version of the relevant h x (t) in the upcom- 
ing subsection. However, we do display all the formulae 
required to construct various types of circular templates 
that are employed in the present study. 



A. GW phasing for compact binaries in 
inspiralling PN-accurate eccentric orbits 

The dominant quadrupolar contribution to the plus 
polarization, denoted by h x \ n (t), associated with a com- 
pact binary consisting of individual masses mi and ni2 
at a radial distance R' [l(j reads 



h x (r(t),<l>(t),f(t),j>(t))\ 



2Gmr]C 
c 4 R' 



Gm 

r{tf I sin2<£(i) 



r{tfj>{tf 



2f(t)r(t)<j)(t) cos20(i) 

'(1) 

where the symmetric mass ratio ?/ = mi mi / m 2 and 
m = m\ + m-2, respectively, while C = cosz, i being 
the inclination of the orbital plane with respect to the 
plane of the sky. The orbital separation, phase and their 
time derivatives are denoted by r(t), <fi(t),r(t) = dr(t)/dt 
and 4>(t) = dcf)(t)/dt, respectively. 

The phasing formalism, detailed in Ref. [l(j, provides 
an efficient way of implementing the temporal evolution 
for the dynamical variables, r(t), 4>(i), rjt), 4>(t), appear- 
ing in Eq. flj. Following Ref. we note that in the 
case of compact binaries in inspiralling eccentric orbits, 
however small may be the eccentricity, the time evolution 
of the dynamical variables explicitly depends on both the 
conservative and the reactive contribution to the orbital 
dynamics. Due to the existence of a Keplerian type para- 
metric solution to the 3PN-accurate conservative orbital 
dynamics, detailed in Ref. [l3T |. the conservative contri- 
butions to |r(i), 4>(t), r(t), </>(£) j can be written as 



r(t)=r (u(l),S,J), 



(2a) 
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r(t)=f (u(l),£,J), (2b) 
4>(t) = £, J) + W(u(l),£, J) , (2c) 
i>(t)=<j>(u(l),£,J), (2d) 

where u and I are the eccentric and mean anoma- 
lies associated with the 3PN-accurate generalized quasi- 
Keplerian parametrization, respectively, while £ and J 
stand for the orbital energy and angular momentum, re- 
spectively. The angular variable <fi is split into two parts 
to make sure that there exists a part that depends lin- 
early on I. Recall that such a split allows one to explicate 
easily the effect of the periastron advance, while con- 
structing the frequency spectrum associated with h x (t) 
relevant for c omp act binaries in PN-accurate eccentric or- 
bits [see Ref. [16| for details]. The angular type variables 
A and I can be expressed as 

I = n(t - t ) + ci = l(u;£,J, q) , (3a) 
A = (1 + k)n(t - t ) + c A = X(t - t ; £, J, c x ) , (3b) 

where the constants to, q and c\ refer to some initial 
instant and the values of / and A at t = to. The sym- 
bols n and k represent the mean-motion and the mea- 
sure of advance of the periastron in the time interval 
2 7r/n, respectively. It turns out that the PN-accurate 
expressions for n and k in terms of £ and J are the 
only two gauge-invariant quantities present in the PN- 
accurate Keplerian-type parametric solution [l3|, [I?) • If 
one ignores the reactive contributions to the orbital dy- 
namics, the explicit temporal evolutions for the dynami- 
cal variables, r{t),(j){t),f(t), (j)(t), are obtained by solving 
the 3PN-accurate Kepler equation (KE) that provides 
the link between I and u. The 3PN-accurate KE can be 
symbolically expressed as 

I =u - e t sinw + h,3(u, £, J) , (4) 

where h,3(u, £,J) stands for the 2PN and 3PN correc- 
tions to the usual Newtonian KE I = u — et sin it. With 
the help of Eq. , we observe that the PN-accurate con- 
servative dynamics, for example, specified by the 3PN- 
accurate relative acceleration A 1 in harmonic gauge [Hj], 
have four constants of integration, namely £ , J , ci and 
ca- 

To describe the dynamics of compact binaries in inspi- 
ralling eccentric orbits, it is essential to include the re- 
active contributions to A % that first enters at the 2.5PN 
(absolute) order. However, it is imperative to invoke the 
energy and angular momentum balance argument to im- 
plement the higher order PN-accurate reactive evolution 
of an eccentric orbit. This is because the reactive con- 
tributions to A % are available only to the 1PN relative 
(or 3.5PN absolute) order [l^]. Recall that the balance 
argument equates the orbital averaged higher order PN- 
accurate energy and angular momentum fluxes in GWs 
to the time derivatives of PN-accurate orbital energy and 
angular momentum, respectively. Invoking the balance 
argument and 2PN-accurate far-zone energy and angu- 
lar momentum fluxes, available in Refs. it 



is possible to incorporate the 2PN-accurate secular RR 
effects on PN-accurate eccentric orbits in the following 
way [see Refs. [HI EH for details]. To describe the inspi- 
ral dynamics of compact binaries, Ref. [Io| employed an 
improved 'method of variation of constants' to include 
the RR effects onto the conservative orbital dynamics. 
The idea is to allow the dynamical variables, appearing 
in Eq. ([T]), have the same functional forms as given by 
Eqs. ([2]) and ([3]) even when the dynamics is not conserva- 
tive, but reactive. However, one lets the associated con- 
stants of integration, namely £ , J, ci and c\ to vary with 
time. The equations that govern the temporal evolutions 
of these quantities are given by Eqs. (35) in Ref. [l(| and, 
as expected, depend on the reactive contributions to A 1 . 
Further, the temporal variations of these four constants 
of integration can be modeled to consist of a slow drift 
and fast oscillations, which symbolically reads 

c a (l) = c a (l) + c a (l) , (5) 

where a stands for one of the four quantities £, J , ci and 
c\. In the above equation, c a (l) denotes the slow drift, 
which accumulates over the RR time scale to induce large 
changes in c a (l). The fast oscillations in c a (l) are de- 
noted by c a {l), and the explicit expressions for dc a /dt 
arc only available to the relative 1PN order [l^]. To 
the reactive PN order we are interested, it turns out 
that dci/dt = dcx/dt = 0, while d£/dt and dj/dt are 
provided by the far-zone energy and angular momentum 
fluxes, respectively [loj . Further, it was demonstrated in 
Refs. [10|, LL5| that the rapidly oscillating parts have sub- 
stantially smaller amplitudes. Therefore, for the present 
study, we ignore the RR induced c a (t) contributions to 
the orbital dynamics. 

The strategy to obtain /i x L(i) associated with com- 
pact binaries moving in 3PN-accurate eccentric orbits 
that are perturbed by the relative 2PN-accurate RR ef- 
fects is the following. With the help of the 3PN-accurate 
generalized quasi-Keplcrian parametrization in harmonic 
coordinates, available in Ref. [l3[, and Refs. [l(| OH, 
we compute the 3PN-accurate expressions for r, W, r, <p 
in terms of £, J and u. This is supplemented by the 
3PN-accurate expression for A in terms of l,£,J. We 
employed a modified version of Mikkola's method [22| . 
detailed in Ref. [1] , to solve accurately and efficiently the 
3PN-accurate KE and obtain the temporal evolution for 
r(£), W(t),r(t),<j)(t) under the 3PN-accurate conservative 
orbital dynamics. Afterwards, we numerically solve the 
2PN-accurate differential equations for £ and J and im- 
pose the resulting reactive evolutions in £ and J onto 
the earlier computed 3PN-accurate conservative orbital 
dynamics. The resulting r(t),(j)(t),r(t), (j>(t) leads to, via 
Eq. (TTJ, /i x | (t) associated with non-spinning compact 
binaries inspiralling along 3PN-accurate eccentric orbits 
under the action of 2PN-accurate reactive dynamics. 

Until now, we have employed the orbital energy and 
angular momentum to describe the eccentric inspiral. 
However, a measure of the non-circularity is better de- 
scribed in terms of any of the eccentricity parameters, 



4 



appearing in the PN-accurate Keplerian-type parametric 
solution |13| . than in terms of J . Therefore, from here 
onwards, we employ the orbital binding energy £ and the 
time eccentricity parameter e*, associated with the PN- 
accurate KE, to specify our PN-accurate eccentric bina- 
ries. With the help of Ref. [Tij . it is fairly straightfor- 
ward to obtain 3PN-accurate expressions for r, r, 4>, W, A 
and 3PN-accurate KE in terms of £ and e t . For the sake 
of brevity, we provide explicitly the above expressions 
to 1PN order only. However, we do provide the 3PN- 



accurate relation connecting n to £ required to obtain 
the lengthy 2PN and 3PN corrections to r, f, <fi, W, A and 
l(u) in terms of £ and e t from Eqs. (23)-(27) in Ref. [U]. 

It should be noted that the dynamical variables en- 
ter h x \~(t) in certain combinations like S^, ^ and -. 
Therefore, we list below the following dimensionless dy- 
namical variables required to describe the dynamical evo- 
lution of /i x |„(/j) in a partially symbolic manner. 



-2~ (£,e t ,u) 
c z r 



- (l + e J" [X(9 - 5v) + 6r? - 16] + e(-) + e (•■) 
X { 4x 



r C 7 e* smu 

-U,e t ,u) = < 1 

c X 



3 9 
8 ~ 8 V 



e+e 2 (...) + e 3 (-) 



(6a) 
(6b) 



r«, e t , „) X ^, e t , «) _ ^ 2 ^ f j i _ + ^ _ + ^ + ^ 



X 



4x 



7T ^-2y [ - 23 + 977 + e 2 (15 - r?) + ^ (^32 -877 + e?(-32 + 8r?)) + £ 2 (...) + £ 3 (■■) ; 



et, it, /) = A(Z, £, e 4 , ) + W(f , e t ,u), where 



G m 



^^(dbfy 



+ e 2 (...) + C 3 (»)K*-*o) 



e t ,u)|ipiv = (v — u) + e t sinu + £ 



(1-e?) 



{(u-u) + et sinit}+£ 2 (...)+£ 3 (..) 



(6c) 
(6d) 

(6e) 
(6f) 



where £ = -2£//jc 2 , /j, being the reduced mass of the where = (1 - Jl - efy/e^. Using the PN accurate 
binary: /1 = 7? to and x — (1 — e t cosu). To evaluate 
W(£, et, it), we employ the following expression for v — u 



v — u = 2 tan 



1 Sill u 



(7) 



relation connecting e^ to e t , extractable from Ref. [1 31 ] . 
it is fairly straightforward to express /3<^> in terms of et 
and £. To 1PN order, /3<^ explicitely reads 



1 + C 



(4 - g) 



+ £ 2 (-)+£ 3 (--) 



(8) 



r 



The explicit expressions for the 2PN and 3PN con- 
tributions to Gm/c 2 r,r/c, X,W,r <j)/c and fy, symboli- 
cally noted as '+£ 2 (...) + £ 3 (..)' in their respective 1PN- 



accurate expressions, are obtainable in a straightforward 
manner using Eqs. (23)- (26) in Ref. [l5[ and the following 
3PN-accurate expression connecting n to £ and et p~3| 



3/2 



Gm 









H 


15 — 7? 


s 128 



(555 + 307?+ 1 It? 2 ) 



(960 - 384r?) 



3072(1 - e 2 ) 3 / 2 



16 (-2340 + {10244 - 123tt 2 }7? - 792r? 2 - 36 e 2 {255 - 159r? + 34r? 2 }) 
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(1 



s 2 ) 3/2 (45 (-653 - 111/7 - lr] 2 + 3r) 3 )) 



The explicit conservative temporal evolution for the 
above listed orbital variables are obtained by solving the 



J 



(9) 



3PN- accurate KE connecting u to I. The right hand side 
of the PN-accurate KE reads 



l(£,u, et) = u — etsinu + ^ 



•15 
, 2 



3r]) (v — u) ^ e t (15 — 77)77 sin(u) 



r 



+ e 3 (-) 



(10) 



where the v — u term appearing at the 2PN order on 
the right hand side of Eq. (fl"0|) is evaluated using 1PN- 
accurate expression for fy. As noted earlier, for an ac- 
curate and efficient determination of u(l), we employed 
a modified version of Mikkola's method for solving the 
classical KE J22j(see Ref. @ for its implementation at 
the 2PN order). With the help of PN-accurate u(l), it is 
fairly easy to obtain h x n W associated with compact bi- 
naries evolving under 3PN-accurate conservative orbital 
dynamics. 

Let us now impose the effects of 2PN-accurate RR onto 
the above detailed conservative orbital dynamics. Due 
to the neglect of the small amplitude and fast oscilla- 
tory contributions to £, J , ci and c\, the reactive orbital 
evolution is entirely given by 2PN-accurate orbital aver- 
aged expressions for the energy and angular momentum 
fluxes, and these PN-accurate expressions in harmonic 



J 



gauge can be extracted from Refs. [19J, [20|, |2l| . However, 
in this study, the PN-accurate eccentric orbits are spec- 
ified in terms of £ and e t , and hence require the 2PN- 
accurate differential equations for £ and et, written in 
terms of £ and e*. Due to the definition of £, the differen- 
tial equation governing the reactive evolution of £ follows 
directly from the 2PN-accurate far-zone energy flux. The 
differential equation describing the secular 2PN-accurate 
reactive evolution of et, requires i) the 2PN accurate ex- 
pression for e 2 in terms of £ and J in harmonic gauge, 
available in Ref. [Hj], ii) the 2PN-accurate expressions 
for the orbital averaged far-zone energy and angular mo- 
mentum fluxes and finally iii) the energy and angular mo- 
mentum balance argument. The resulting 2PN- accurate 
expressions for d^/dt and det/dt, describing the secular 
evolutions of £ and e< , is given by 



d£ = 
dt 

de t _ 
dt 

c 



~1PN 



A1.5PN 



i'2PN 



15 Gm 

£ 4 e t 77 c 3 



JV _|_ £1PN 



e t 



il.&PN 



15 Gm 

{96 + 292 e 2 + 37e^} 



£ 2PAr }, 

e 2PN ) , where 



(I-?) 772 
£ 

2N9/2 



56 (1 - e t 

e /2 Tr 



{208 - 13440?7 + 8[22340 - 1960777]e 2 + 42[5966 - 3459?7]e 4 + [19487 - 880677]^} 



90 



{36864 + 448320e t 2 + 2061840e 4 + 6204647e?} 



e 



— {32 [253937 - 162585?; + 4536077 2 ] + 16 [2007326 - 6275466?7 + 350954177 2 ] e 



6048(^1 - e 2 ) 1 

- 12 [22179502 - 3510340277 + 13897611?7 2 ] e 4 + 18 [12638659 - 12297567?? + 3878616?7 2 



(lla) 
(lib) 
(11c) 

(lid) 
(lie) 



+ 9 [1374197 - 916728?7 + 254856?7 2 ] ef - 9072 [96 - 1060e 2 - 1863e 4 - 148ef] ^1 - e 2 (5 - 2?/)} , (llf) 
N _ (l21e 2 + 304) 



(Hg) 
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• 1PN 



l.bPN 



2PN 



56(yr^fr 
e /2 7T 



{75667^ + 344784e 2 - 980 (34^ + 225e 2 + 72) 77 - 28536} 



96 



189120 + 1042992 e t 2 + 3061465 e t A 



(llh) 

(in) 



i 2 



: 17914878- 12062097 jy+ 3288852 yy 2 ) e t 6 



2016(v/T^e|) 9 

(128983907 - 139206105 r, + 45633420 rf) e t 4 + (46339524 ?/ 2 - 94111416 r, + 28071312) e t 2 



+ 7116856 - 6772968 77 + 3566304 if + 3024 (511 e t 4 + 2809 e t 2 + 80) (5 - 2 jj) J 1 - e 



(Hj) 



The contributions entering at the relative 1.5PN order 
are due to the dominant order GW tails that are usually 
expressed in terms of certain infinite sums involving the 
Bessel functions J p (pe t ) and its derivative J' (pet) fl9l . 
[20| . In the above equations, we have expanded those 
infinite sums to 0(ef). This is justified as we are only 
interested in eccentric binaries having initial e t < 0.2 in 
the present study. Let us remind the reader that it is 
due to the neglect of the tilde contributions to £ 7 e t ,ci 
and c\ that the reactive evolution is fully prescribed by 
the above differential equations for £(i) and e*. And, the 
above two equations represent the secular variations in 
£ and e t , namely £(t) and e~t(t), due to the emission of 
GWs. 

We are now in a position to implement numerically 
h x \ n (t) emanating from non-spinning compact binaries 
inspiralling along 3PN-accurate eccentric orbits under 
the action of 2PN-accurate reactive dynamics. It is ob- 
vious that we need to specify the initial and final val- 



ues for £ so that the dominant GW spectral compo- 
nent of /ix|q(*) evolves in the frequency window 40Hz 

- (6 3/2 7rGm/c 3 ) _1 Hz, relevant for the initial LIGO. 
We compute the bounding values for £ with the help 
of Rcf. [16| that demonstrated the fact that the domi- 
nant GW spectral component of a mildly eccentric com- 
pact binary, having PN accurate orbital motion, appears 
at i < 4r >= l£±£ where < 4r > stands for the 

tt at 7r ' at 

PN-accurate orbital average of d<j)/dt. The PN-accurate 
expression for < > can be computed using the proce- 
dure, detailed in Ref. [23[, involving the Laplace second 
integrals for the Legendre polynomials [see Eqs. (4.16)- 
(4.19) in Ref. (23|] . For the present paper, we use the fol- 
lowing 3PN-accurate expression for the orbital averaged 
orbital angular velocity < -± > to obtain the bounding 
values for £ for a compact binary specified by certain m, 77 
and et values: 



dt 



n(l + k) 

r.3 ( 



3 C 



15 



128(1 



(1851 - 786 77 + 11 ?7 2 ) - 2 e 2 (-861 + 486 77 + 11 t? 2 
3 (76965 + (-167089 + 3936tt 2 ) 77 



+4(555 + 30 ?7 + 11 V 2 ) + 192(1 - e 2 ) 3 / 2 (-5 + 277) 

+5439 ?7 2 + 45 77 s ) + 3 e 2 (228249 + (-222781 + 984tt 2 ) 77 + 40491 t? 2 - 135 77 s ) + 9 ^(25301 - 17201?7 + 3471?7 2 
+45?7 3 ) - 45e^(-653 - 111?7-7t7 2 

(12) 



3»7 3 ) + V 1 - e 2 (l6(-6660 - 41(-292 + 3tt 2 )?7 - 792 v 2 
+e 2 (6120 + (-9704 + 1237r 2 >7 - 432?7 2 ) + 36et (255 - 159^ + 34?7 2 ))) 



With these inputs, it is fairly straightforward to construct 
h x \ n (t) for inspiralling compact binaries, having some 
tiny orbital eccentricity, emitting GWs in the initial- 
LIGO frequency window. 



Finally, let us explain why we employed £, along with 
et, to represent the PN-accurate eccentric orbit, which 
guided us to use £ as our PN expansion parameter. In the 
literature, the GW phasing for eccentric binaries are per- 
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formed using either (G™«) (2/3) [n| |i| or (G™h) (2/3) 
jMEl, where to stands for the above listed 3PN-accurate 

< ^ >, as the PN expansion parameters, along with et- 
The choice of n = 2n/T ri where the radial period T r 
being the time of return to the pcriastron, as a PN ex- 
pansion parameter may be visually problematic. This 
is because it is rather difficult to pinpoint the pcrias- 
tron position for circular binaries, making it difficult to 
imagine the meaning of T r in the case of circular inspi- 
ral. We speculate that this may be the one of the rea- 
sons for Ref. m to introduce to =< ^ >= n (1 + k) 
as their PN expansion parameter. However, it is not 
clear why one should employ the orbital averaged an- 
gular velocity to perform GW phasing for eccentric bi- 
naries. The reason stated in Ref. [24| for introducing 

< ^ > as the PN-expansion parameter is that it gener- 
alizes the to = d(j)/dt parameter used in the construction 
of the usual PN-accurate inspiral templates, available in 
Ref. In other words, the PN accurate expression 
for d4>/dt, given in terms of £, et,u, when re-expressed 
in terms to =< 4^ >,e t ,u using Eq. (fT2")l will lead to 
d<fi/dt = to in the circular limit [check Eqs. (A.fO)-(A.14) 
in Ref. [25[ f° r the actual demonstration of the above 
statement]. It is important to realize that the use of 
d<p/dt = to as the PN expansion parameter for the cir- 
cular inspiral is usually justified by the argument that it 
defines the instantaneous orbital angular frequency. In 
the case of the heavily employed restricted PN-accurate 
inspiral templates [J], the instantaneous GW frequency 
/gw is related to the instantaneous orbital angular fre- 
quency by /gw = to/ir. However, in the case of GW 
phasing for eccentric binaries, neither n nor lo =< ^ > 
are sufficient to define the instantaneous orbital angu- 
lar frequency. Further, the associated GW spectrum can 
not be expressed solely in terms of n or < ^ > Q. 
These considerations prompted us to employ neither n 
nor lu =< ^ > as our PN expansion parameter for per- 
forming the GW phasing for eccentric binaries. We would 
like to stress that, while considering the PN-accurate con- 
servative orbital dynamics, our PN expansion parame- 
ter requires no orbital averaging or any prior knowledge 
about the turning points of the associated eccentric or- 
bit. Further, £ provides a measure of the instantaneous 
orbital energy of an eccentric orbit. In the case of cir- 
cular inspiral, it is also straightforward to identify £ as 
a measure for the instantaneous orbital energy. Interest- 
ingly, for Newtonian circular orbits, our PN expansion 
parameter £ is indeed a measure of ^ like ( G ™" ) 2 / 3 
and ( G "" ) 2 / 3 variables available in the literature. 

It was argued in Ref. [T(| that the GW phasing for 
eccentric binaries that employ n as the PN expansion 
parameter is not going to be accurate near the last sta- 
ble orbit (LSO). This is one of the reasons for termi- 
natin g th e eccentric orbital evolution around j 2 ~ 48 in 
Refs. [19, 0, wherc 3 = C 'J/G 

mfi. Recall that for a 
test particle around a Schwarzschild black hole, the last 
stable circular orbit is also given by j 2 = 12 [17|. A re- 



cent comparison between PN and Numerical Relativity 
(NR) based GW phase evolutions for equal-mass eccen- 
tric binaries, reported in Ref. |25|l. also supported the 
above argument. However, Ref. [251 ] demonstrated that 
the PN-accurate GW phasing for eccentric binaries, hav- 
ing 3PN-accurate conservative and 2PN accurate reactive 
dynamics, that employ to =< > as the PN expansion 
parameter arc fairly accurate even near the LSO. This 
observation can be linked to the fact that the GW phas- 
ing for eccentric binaries in terms of to =< > and e t , 
in the limit e t — > provides the 2PN-accurate TaylorT4 
approximant. And, Refs. [H, H^l recently demonstrated 
that the GW phase evolution under this TaylorT4 ap- 
proximant is fairly close to its NR counterpart for equal 
mass binary black holes having around 10-15 orbital cy- 
cles close to their LSOs. We would like to point out that 
in the limit et — > 0, our present prescription for doing 
GW phasing for eccentric binaries leads to the TaylorEt 
approximant that employ the 3PN accurate expression 
for d<p/dt and the 2PN accurate expression for d£/dt [ill ]. 
With the help of Fig. 4 in Ref. (26|, we can state that the 
the accumulated GW phase difference between the above 
approximant and NR, during the late inspiral stage of an 
equal-mass binary, is comparable to what is provided the 
TaylorT4 approximant at 2.5PN order. The above phase 
difference turned out to be substantially smaller than 
a similar accumulated GW phase difference originating 
from a 2PN approximant that employs ( G "" ) 2 ^ 3 as the 
PN expansion parameter [Til ]. Therefore, it is reasonable 
to expect that during the late inspiral stage, our present 
prescription for doing eccentric GW phasing could be 
closer to NR than what is presented in Refs. [TO, 111]- We 
are initiating an effort to make such a comparison with 
the the help of yet-to-be available fully 3PN-accurate RR 
effects for the eccentric inspiral, similar to what is re- 
ported in Ref. [HI- 

In what follows, we briefly describe how we construct 
various PN-accurate circular inspiral templates that in- 
corporate RR effects to the relative 2PN accuracy. 



B. Various PN-accurate circular inspiral templates 
that incorporate RR effects at the relative 2PN 
order 

In this subsection, we briefly describe how to construct 
the earlier mentioned PN-accurate templates relevant for 
the circular inspiral. The crucial PN-accurate quantities 
that are required to construct these circular inspiral tem- 
plates are the 3PN-accurate orbital binding energy £(x ) 
[Hm and the 2PN-accurate GW luminosity C{x) J28[, 
both of which are usually expressed as PN series in the 
gauge invariant quantity x = (Gtooj/c 3 ) 2 / 3 , where to is 
the instantaneous orbital angular frequency of the asso- 
ciated circular binary. The explicit expressions for £ (x) 
and C(x) relevant for the present study are 
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The time-domain x-based circular templates are also 
constructed under the restricted PN-accurate prescrip- 
tion and this leads to 



AGmrjC 
c 2 R' 



x(t) sm2 (j)(t) , 



(14) 



Following Refs. [1, Q, one can construct different tem- 
plate families that provide different prescriptions for the 
evolution of x(t) and <f>(t) appearing in the above equa- 
tion. The 2PN order adiabatic TaylorTl approximant of 
Ref. [J] is given by 



dtj>(t) 

dt 
dx{t) 

dt 



w(t) 



,.3/2 



Gm 
C{x) 



(d£ 2 p^/dx) ' 



(15a) 
(15b) 



where ^2PN stands for the 2PN-accurate orbital binding 
energy, obtainable by dropping x 3 terms in Eq. (|13a[) for 



£{x). It should be noted that in this approximant, avail- 
able in LAL, one truncates the PN-accurate expressions 
for £ (x) and C(x) at the same relative PN order. We can 
also create the complete adiabatic TaylorTl approximant 
in the following way Q 



dcj>(t) 

dt 
dx{t) 

dt 



w(t) = 



,.3/2 



Gi 



(d£/dx) 



(16a) 
(16b) 



where £{x) is indeed 3PN-accurate and is given by 
Eq. (fT3aj) . 

The Taylor T4 approximant, introduced in Ref. (l2l |. is 
obtained by Taylor expanding in x the right-hand-side of 
either Eq. (|15bj) or Eq. (|16bj) to the 2PN order in x. The 
2PN-accurate TaylorT4 is defined by 



^="(*) = t4* 3/2 ' ( 17a ) 

dt Gm 

dx{t) c 3 64r/ .( /743 11 \ A , /2 /34103 13661 59 2 \ 2 \ 

— = -x 5 il-{ 1 T7 } x + 4irx A/2 + [ 1 r?H n 2 ) x 2 > . (I7b) 

dt Gm 5 \ \336 4 ') \l8lte 2016 ' 18 ' J / v ; 



While searching for inspiral GWs, the LSC employs (or 
has employed even in the recent past) the Fourier domain 
2PN-accurate TaylorF2 approximants ||. This approxi- 
mant can be constructed analytically from the above Tay- 
lor T4 a ppr oximant using the stationary phase approxi- 
mation [29| making it computationally cheaper than the 
time-domain TaylorT4 approximant. However, the fo- 
cus of the present study is not the computational costs 
associated with the implementation of various template 



families, therefore, we employed the 2PN-accurate time- 
domain TaylorT4 approximant rather than its Fourier 
domain counterpart. We observe that among the three 
approximants, the LSC Algorithm Library (LAL)[3(| 
strictly provides only one of them, namely, the adiabatic 
TaylorTl approximant at the 2PN order. 

Finally, we construct the 2PN version of the TaylorEt 
using Ref. [ll| as 



h{t) 

dcj)(t) 

; = 

dt 



AGmrjC 
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m sin 2 0(t) 
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The above approximant obviously provides the circu- 
lar limit of GW phasing for eccentric inspiral, detailed in 
the previous subsection, and hence employ £ as the PN 
expansion parameter. A number of features of this ap- 
proximant are worth mentioning. The TaylorEt approxi- 
mant is specified by two separate PN-accurate differential 
equations and in contrast, the various x-bascd approxi- 
mants arc governed by just one PN-accurate differential 
equation. In a given GW frequency band, the limits of 
the integration for various x-based approximants, irre- 
spective of whether one is interested in the Newtonian 
or 3.5PN order approximant, are the same and are in- 
dependent of the symmetric mass ratio rj. However, the 
limits of the integration for the TaylorEt approximant 
do depend on the underlying 3PN-accurate conservative 
dynamics and therefore depend on rj. 

The plots of H x (t) representing h x \ Ci (t) scaled by 

GmrjC / c 2 R' evolving under various prescriptions de- 
tailed in this section are displayed in Fig. [TJ These plots, 
displaying the evolving H x (t) in the initial-LIGO fre- 
quency window, represent the canonical neutron star- 
black hole binary. The plots are for i) an eccentric binary 
inspiralling with an initial et = 10 -3 , ii) the TaylorEt ap- 
proximant, iii) the adiabatic TaylorTl approximant and 
iv) the TaylorT4 approximant. The similar duration for 
the eccentric and the circular TaylorEt H x (t) should be 
noted. 

Let us now move on to explore the ability of various 
circular templates to capture GWs from inspiralling com- 
pact binaries having tiny orbital eccentricities. 



III. PERFORMANCES OF VARIOUS 
CIRCULAR TEMPLATES AGAINST GW 
SIGNALS FROM COMPACT BINARIES IN 
INSPIRALLING ECCENTRIC ORBITS 

In this section, we explore how effectual and faith- 
ful are the various 2PN order circular templates, intro- 
duced in the previous section, with respect to GW sig- 
nals from compact binaries in inspiralling eccentric orbits 
that should be interesting to the initial-LIGO data ana- 
lysts. This is achieved by computing the fully maximized 
overlaps involving the eccentric h x (t), treated to be the 
fiducial inspiral GW signal, and the four types of circu- 
lar inspiral templates, given by Eqs. (fl5l) . (fTBj) . (fT?) and 
([TB| . We extensively employ Ref. [l4[ to compute the 
fully maximized overlaps with the help of their minimax 
overlap. In our numerical experiments, the maximization 
over the time-lags are straightforward as we invoke the 
fast-Fouricr-transform routine of the Numerical Recipes 



j3l| to compute the Fourier-domain versions of our time- 
domain GW signals and templates. Further, the maxi- 
mization over m and ?y, characterizing the various circular 
templates are performed using the amoeba routine 31| . In 
the present context, a given circular template family is 
considered to be effectual and faithful if some members of 
the family can provide fully maximized overlaps, hereon 
the fitting factors (FFs), that are closer to unity and have 
smaller biases in estimating the parameters representing 
our fiducial GW signal. In other words, the values of 
m and rj characterizing the circular template having the 
largest FF provide a measure of the faithfulness of that 
particular template family. It turns out that the drop in 
event rate is proportional to (1 — FF 3 ) [32j and therefore, 
it is indeed desirable to have circular template families 
having FF values that are > 0.97 with respect to fidu- 
cial GW signals from compact binaries having tiny or- 
bital eccentricities. This is mainly because the expected 
astrophysical GW sources will have tiny, but non-zero, 
orbital eccentricities. 

In Table [IJ we show the FFs involving our fiducial 
eccentric GW signals and the 2PN order Taylor T4 cir- 
cular templates, relevant for the initial LIGO. We re- 
strict our attention to the usual canonical binaries rep- 
resenting neutron star-neutron star, black hole-neutron 
star and black hole-black hole binaries. In our FF 
computations, we evolve both our fiducial GW signals 
and circular templates in the frequency range 40 Hz to 
(6 3 / 2 7r G m/c 3 ) Hz. The numbers clearly indicate that 
the 2PN TaylorT4 circular templates arc neither effec- 
tual nor faithful in extracting our fiducial GW signals 
from inspiralling binaries, having low orbital eccentrici- 
ties like e t ~ 0.001. The numbers associated with the 
e t = 10~ 3 case are representative of what to expect for 
still lower values for the initial et- We justify the above 
statement by repeating the FF computation employing 
the circular TaylorEt approximant to model the fiducial 
GW signal. And, the resulting numbers are listed in the 
first row of Tabic fl] 

Next in Table [Til wc show the results of our FF com- 
putation involving the adiabatic and complete adiabatic 
TaylorTl circular templates that require 2PN accurate 
RR effects. The complete adiabatic TaylorTl templates 
are invoked so that we can include the effects of 3PN- 
accurate conservative orbital dynamics into the usual x- 
based circular templates that require the 2PN-accurate 
reactive dynamics, given by Eqs. Q13p . The circular Tay- 
lorEt approximant is again invoked to model our fiducial 
GW signal mainly to point out that the FFs are not im- 
proved for binaries having initial et < 0.001 . 

The performances of the TaylorEt circular templates 
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against our fiducial GW signals from eccentric binaries 
are compared in Table IIIH The numbers, not surpris- 
ingly, demonstrate that the circular inspiral templates 
based on the TaylorEt approximant is both effectual and 
faithful in capturing GWs from compact binaries having 
tiny orbital eccentricities. The listed biases in the total 
masses of the canonical binaries are consistent with the 
fact that both eccentricity and higher values for the total 
mass are capable of decreasing the number of GW cycles 
in a given GW frequency window. 

The numbers displayed in Tables Q] and [IT] clearly indi- 
cate that GW signals from inspiralling compact binaries 
having tiny orbital eccentricities, constructed by adapt- 
ing the phasing formalism, detailed in Ref. [n| EH an( i 
summarized in Sec. [Til can not be captured by the usual 
x-based adiabatic and complete adiabatic TaylorTl and 
adiabatic TaylorT4 circular templates that incorporate 
the 2PN-accurate reactive dynamics. This is should be 
of definite concern for the GW community because the 
expected GWs from astrophysical compact binaries will 
have tiny orbital eccentricities and the GW data ana- 
lysts belonging to the LSC employed, even in the recent 
past, the traditional x-based 2PN order circular inspi- 
ral templates to search for GWs from non-spinning com- 
pact binaries Q . The numbers displayed in Tables [J and 
HI! evince that the dominant reason for the above unex- 
pected conclusion is the inability of the x-based circular 
inspiral templates, based on the two TaylorTl and T4 
approximants, to capture GW signals modeled after the 
TaylorEt approximant, in an effectual and faithful man- 
ner. 

From a purely mathematical point of view, it may be 
argued that the circular TaylorEt and x-bascd approxi- 
mants should be identical as they differ only in the use 
of the PN-expansion parameter. This implies that the 
observed differences in GW phase evolutions under var- 
ious PN approximants may be attributed to the uncon- 
trolled higher order terms. However, the present anal- 
ysis and the one detailed in Ref. [33[ indicate that for 
the data analysis considerations the inspiral templates 
belonging to the TaylorEt and the x-bascd approximants 
should be treated to be different. From a physical point 
of view, these approximants are also different due to the 
following reasons. It should be noted that due to the 
use of the standard energy balance argument, various 
circular inspiral templates model GWs from compact bi- 
naries evolving along a sequence of circular orbits under 
the action of gravitational radiation reaction. For the x- 
based templates, the above circular orbits can be treated 
to be Newtonian-accurate due to the use of Eq. Q15ap 
for d(f>/dt. However, the Taylor Et approximant requires 
the 3PN-accurate expression for d<p/dt and therefore the 
sequence circular orbits can be considered to be 3PN- 
accurate. The differences in the above mentioned circular 
orbits imply that the TaylorEt approximant is specified 
by two separate PN- accurate differential equations, while 
the TaylorTl and TaylorT4 approximants are governed 
by just one PN-accurate differential equation. Further 



and as noted earlier, for a given GW frequency win- 
dow, the limits of the integration for the x-based ap- 
proximants, irrespective of whether one is interested in 
Newtonian or 3.5PN accurate templates, are the same 
and are independent of r\. However, the limits of the 
integration for the TaylorEt approximant in the present 
study crucially depend on the underlying 3PN-accurate 
conservative dynamics and do depend on r\. 

It may be argued that if one constructs eccentric GW 
signals using the generalized x parameter, introduced in 
an ad- hoc manner in Refs. (24l.[25j. and et, the 2PN order 
TaylorT4 approximant (and hence LSC's 2PN-accurate 
TaylorF2 approximant) will have no difficulty in cap- 
turing such eccentric GW signals. This is because the 
TaylorT4 approximant indeed provides the circular limit 
of GW phasing for eccentric binaries that employ the 
generalized x parameter and et. Let us (again) recall 
that the use of u) = d(f>/dt as a PN expansion param- 
eter for the circular inspiral is justified due to the fact 
that it provides the instantaneous orbital angular fre- 
quency. Further, the emitted instantaneous GW angular 
frequency is related to its orbital counterpart by a fac- 
tor of two in the case the restricted PN-accurate circular 
templates. However, as noted in the previous section and 
clearly visible in Eqs. (A10)-(A14) in Ref. [HI, the PN 
expansion parameter uj =< d(f>/dt > present in the gen- 
eralized x parameter will not provide the instantaneous 
orbital angular frequency in the case of compact binaries 
in inspiralling eccentric orbits. Further, the GW spec- 
trum associated with compact binaries in PN-accurate 
eccentric orbits is not expressible as integer multiples 
of < d<p/dt >= nil + k) Q. Therefore, it is not clear 
why one should employ u) =< d<fr/dt > as the PN ex- 
pansion parameter to perform GW phasing for eccentric 
binaries, when there exists a more primary and gauge- 
invariant variable like £ that is related to the binary's 
instantaneous orbital energy. 

The practitioners and proponents of the x-based in- 
spiral h x _ + (t), both circular and eccentric, may point 
out that the PN-accurate GW phasing for eccentric bi- 
naries in terms of the generalized x parameter and et are 
closer to what is provided by NR, as reported in Ref. [25| . 
In our opinion, this is most likely due to the ability 
of the circular 2PN-accurate TaylorT4 approximant to 
track fairly closely the GW phase evolution, originating 
from NR, associated with the equal-mass binary black 
hole inspirals that last several orbits near their respec- 
tive last stable orbits [l^, HI]. An inspection of Fig. 4 in 
Ref. (2|| reveals that at the 2PN reactive order the accu- 
mulated GW differences between TaylorEt-TaylorTl and 
TaylorEt- TaylorT4 are ~ 4.4 radians and ~ 2.6 radians, 
respectively. These differences in the accumulated GW 
phase may be attributed to the uncontrolled higher order 
terms present in the various Taylor approximants and are 
partially responsible for our FF results. The statement 
about the partial responsibility is due to our observa- 
tion that both the TaylorTl and T4 approximants are 
fairly effectual in capturing GWs from equal mass circu- 
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lar compact binaries, based on the TaylorEt approximant 
and having to > 30 Mq. However, these x-based approx- 
imants are not effectual (and faithful) in the case of low 
mass binaries. We note that similar results prevail even 
when the RR order is increased to the 3.5PN relative or- 
der while dealing with equal-mass circular inspirals(33T|. 
This is despite the observation that the GW phase evolu- 
tion under the circular TaylorEt approximant at 3.5PN 
order is fairly close to NR (and hence to the TaylorT4 
approximant at 3.5PN order) in comparison with what 
is provided by the TaylorEt approximant at the 2PN re- 
active order |2q |. 

The Fig. 4 in Ref. [26[ indicates that the PN-accurate 
GW phase evolutions prescribed by the circular TaylorEt 
and TaylorT4 approximants at 2.5PN reactive order can 
not be substantially different from each other. Therefore, 
it is natural to expect that these two approximants will 
be effectual and faithful to each other for data analysis 
considerations. The entries listed in Table IIVI support 
somewhat different conclusions, similar to what is de- 
tailed in Ref. [33j ] - In other words, for low to compact 
binaries, the Taylor T4 approximant is not effectual and 
in the case of high to binaries, the x-based approximant 
is effectual at the expense of inviting some what large 
biases in r\. 

The earlier discussions that supported the use of £ as 
an appropriate PN expansion parameter for doing GW 
phasing for eccentric binaries and the monotonic conver- 
gence to the NR based GW phase evolution, exhibited 
by its circular counterpart, suggest that the TaylorEt 
approximant at 3.5PN order is a promising candidate 
to search for inspiral GWs from astrophysical compact 
binaries that are surely going to have tiny orbital eccen- 
tricities. Further, the Fig. 4 in Ref. [26| reveals that the 
accumulated phase disagreement between the 3.5PN or- 
der TaylorEt approximant and 2PN order TaylorT4 ap- 
proximant is ~ 1.00 radians while considering an equal- 
mass binary inspiral that last around nine orbits during 
its late inspiral stage. Therefore, we probed how effec- 
tual and faithful are these two distinct approximants and 
the results are listed in Table [V] We observe that FFs 
that are greater than 0.97 are only possible for binaries 
that have to > 35 Mq and for other stellar mass bina- 
ries, the 2PN order TaylorT4 templates arc fairly inef- 
fectual. It is worthwhile to note that this experiment 
is also prompted by the observation that the LSC em- 
ployed 2PN-accurate TaylorF2 templates that arc identi- 
cal to the 2PN-accurate TaylorT4 approximant for data 
analysis considerations [||. 

It is plausible that the inability of the 2PN-accuratc 
TaylorT4 approximant to capture inspiral GWs modeled 
after the TaylorEt approximant or its eccentric extension 
can have astrophysical implications. Recently, the LSC 
reported that it is rather impossible to have a compact 
binary progenitor for the gamma-ray-burst GRB 070201 
having its masses in the range 1M Q < mi < 3 Mq and 
IMq < to 2 < 40 Mq to be located in the spiral arm of 
the Andromeda galaxy This conclusion was reached 



by analyzing the LIGO data, associated with the sci- 
ence run S5, lasting around 180 seconds around the time 
of the GRB 070201. Interestingly, the LSC employed 
the 2PN-accurate TaylorF2 templates to filter the LIGO 
data while searching for inspiral GWs from non-spinning 
compact binaries in the above mass range [3~i| . It is astro- 
physically quite possible that the above progenitor can- 
didate will have some tiny non-zero residual orbital ec- 
centricity and therefore, we explored the ability of the 
2PN-accurate TaylorT4 approximant to extract the in- 
spiral GWs from compact binaries having residual ini- 
tial e t ~ 10~ 3 . The use of the 2PN-accurate TaylorT4 
approximant is justified as the 2PN-accurate TaylorF2 
templates are only the analytical Fourier-domain version 
of the time-domain TaylorT4 templates. The numbers 
listed in Table IVII indicate that the TaylorT4 approxi- 
mant is effectual with respect to our eccentric GW signal 
only in the high total mass scenario like the case involv- 
ing mi = 3 Mq , to 2 = 40 Mq and in these cases the 
search templates are highly unfaithful. For other mass 
combinations, the TaylorT4 templates are neither effec- 
tual nor faithful. The entries listed in Table IVII force 
us to state that further investigations involving the rele- 
vant LIGO data sets and the various PN-accurate LAL 
templates would be very useful to clarify the real impli- 
cations of our FF results. Finally, we would like to state 
that the above detailed FF experiments were repeated 
several times by employing a wide range of initial val- 
ues for to and 77, required by the amoeba routine, and 
the numbers listed in our tables belong to the runs that 
provided the best possible FFs. 

The results arising from our various numerical exper- 
iments force us to conclude that it is rather desirable 
(and may be urgent) to explore the ability of various 
x-based PN-accurate inspiral templates, prescribed in 
Refs. 0, [13], to capture inspiral GW signals from as- 
trophysical compact binaries that are bound to have tiny 
residual orbital eccentricities. 



IV. CONCLUSIONS 

In this paper, we probed the ability of a number of tra- 
ditional PN-accurate circular inspiral templates, based 
on the adiabatic TaylorTl, completely adiabatic Tay- 
lorTl and adiabatic TaylorT4 approximants, to capture 
inspiral GWs from compact binaries having tiny resid- 
ual orbital eccentricities. Our inspiral GW signals are 
modeled by adapting the phasing formalism of Ref. [T(| 
and therefore originate from compact binaries moving 
in 3PN-accurate eccentric orbits, specified by £ and et, 
that are perturbed by the 2PN- accurate RR effects. We 
demonstrate that the traditional circular templates that 
incorporate 2PN-accurate RR effects are neither effec- 
tual nor faithful in capturing our GW signals having 
tiny orbital eccentricities. However, PN-accurate circu- 
lar templates based on the recently introduced TaylorEt 
approximant are both effectual and faithful with respect 
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to GW signals having small orbital eccentricities like 0.1. 
We explained the physical motivation for employing £, a 
measure of the instantaneous orbital binding energy, to 
describe PN-accurate circular and eccentric orbits. With 
the help of a number of FF experiments, we explored the 
implications of using only the traditional PN-accurate 
inspiral templates to filter the LIGO data. 

We are planning and pursuing a number investigations 
motivated by the physical, numerical relativity and data 
analysis considerations to probe various aspects of the 
TaylorEt approximant, both circular and eccentric. It 
will be desirable to extend the present study by including 
at least the 3PN-accurate RR effects in an accurate and 
efficient manner and which include the c a contributions 
while describing the inspiral dynamics of PN-accurate ec- 
centric binaries, relevant for both the ground-based and 
the proposed space-based GW interferometers. Using 



such a PN-accurate dynamics, it will be interesting to 
reexamine what is pursued in Ref. p5j . while employing 
£ rather than u> =< > to describe the PN-accurate 
orbit. We are in the middle of incorporating the spin ef- 
fects into the present study with the help of inputs that 
generalize what is detailed in Ref. [35| . 
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FIG. 1: A set of plots that display the temporal evolution for H x (t) = /i x IqM/ ( G ^'f ) i originating from a neutron star-black 
hole binary having mi = 10 Mq and 7712 = 1.4 Mq, and in the initial LIGO frequency window. The top left panel depicts an 
eccentric orbital evolution of Sec. [II] having an initial orbital e t = 10 -3 . The other three plots are for the circular templates 
modeled after the TaylorEt approximant, the adiabatic TaylorTl and T4 approximants and these a given by Eqs. (|18|l . (|15[) 
and 1)171) . respectively. The accumulated number of GW cycles in the initial LIGO frequency window, defined by 40Hz- 

^6 3 ' 2 Gni7r/c 3 j Hz, are ~ 346.264,346.257,326.048,333.622 for the eccentric inspiral and circular inspirals defined by the 
TaylorEt, adiabatic TaylorTl and T4 approximants, respectively. 
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TABLE I: The initial LIGO FFs involving the 2PN-accurate circular TaylorT4 approximant for the three types of canonical 
binaries. The fiducial GW signals from compact binaries in inspiralling eccentric orbits are constructed using the phasing 
formalism detailed in Sec. [TTJ The templates that provide the listed fully maximized minimax-overlap (FF) are characterized 
by mt and 774. In the case of double neutron-star binaries, we terminate the the temporal evolution of both the GW signal 
and circular templates when their respective GW frequencies reach 1000Hz. In the et = cases, we employ the TaylorEt 
approximant to model the fiducial GW signals and the associated numbers can be considered to represent what to expect for 
initial et < 10 -3 . The numbers indicate that the 2PN-accurate TaylorT4 templates are neither effectual nor faithful while 
dealing with inspiral GWs from compact binaries having tiny orbital eccentricities. 



TUi/Mq : 7712 /Mq 


1.4 : 1.4 


1.4 : 10.0 


10.0 : 10.0 


e t = FF 
mt 
r\t 


0.822 
2.800 
0.248 


0.939 
6.906 
0.250 


0.957 
19.04 
0.247 


et = 0.001 FF 

mt 
Vt 


0.821 
2.800 
0.248 


0.939 
6.906 
0.250 


0.951 
20.00 
0.224 


e t = 0.01 FF 

m t 
Vt 


0.821 
2.800 
0.248 


0.938 
6.909 
0.250 


0.951 
20.00 
0.225 


e t = 0.10 FF 

mt 
Vt 


0.708 
2.800 
0.250 


0.830 
6.943 
0.250 


0.830 
19.00 
0.250 



TABLE II: Values of the FFs and the associated biases in m and 77 values, while employing the circular inspiral templates based 
on the adiabatic Tl and the complete adiabatic Tl approximants. The complete adiabatic TaylorTl templates are employed 
to incorporate the effects of 3PN-accurate conservative dynamics into the z-based circular templates The conclusions and 
other details are similar to what is observed in Table U 



mi/M e : 7712 /M© 1 1.4 : 1.4 1 1.4 : 10.0 


10.0 : 10.0 


Standard adiabatic TaylorTl templates 


e t = FF 
mt 
Vt 


0.771 
2.785 
0.250 


0.843 
6.845 
0.250 


0.894 
18.55 
0.243 


e t = 0.001 FF 

mt 
Vt 


0.768 
2.800 
0.248 


0.842 
6.846 
0.250 


0.891 
18.57 
0.242 


e t = 0.01 FF 

mt 
Vt 


0.766 
2.800 
0.248 


0.841 
6.834 
0.250 


0.841 
6.837 
0.250 


e t = 0.10 FF 

mt 
Vt 


0.677 
2.800 
0.249 


0.766 
6.876 
0.250 


0.778 
20.01 
0.216 


Complete adiabatic TaylorTl templates 


e t = FF 
mt 

m 


0.761 
2.783 
0.250 


0.800 
6.820 
0.250 


0.878 
18.41 
0.245 


e t = 0.001 FF 

mt 
Vt 


0.761 
2.786 
0.250 


0.800 
6.821 
0.250 


0.878 
18.41 
0.246 


e t = 0.01 FF 

mt 
Vt 


0.758 
2.800 
0.247 


0.800 
6.820 
0.250 


0.878 
18.41 
0.247 


e t = 0.10 FF 

mt 
Vt 


0.673 
2.800 
0.249 


0.754 
6.867 
0.249 


0.838 
18.55 
0.247 
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TABLE III: The results of FF computations that employ the circular templates based on the TaylorEt approximant. It is 
important to note that the temporal evolution for both the eccentric GW signals and the circular templates require the 3PN 
accurate conservative and the 2PN-accurate reactive dynamics The rest of the details are as in Tables U and |TI] The displayed 
values clearly demonstrate, not surprisingly, that the circular TaylorEt templates are both effectual and faithful with respect 
to our inspiral GW signals having tiny orbital eccentricities. 



mi/M® : 7712 /Mq 


1.4 : 1.4 


1.4 : 10.0 


10.0 : 10.0 


e t = 0.001 FF 

mt 
Vt 


0.999 
2.824 
0.246 


I. 000 

II. 50 
0.106 


1.000 
20.17 
0.246 


e t = 0.010 FF 

mt 
Vt 


0.999 
2.823 
0.247 


I. 000 

II. 45 
0.107 


1.000 
20.18 
0.246 


e t = 0.100 FF 

mt 
Vt 


0.804 
2.811 
0.250 


0.921 
11.44 
0.108 


0.956 
20.20 
0.249 



TABLE IV: The initial LIGO FFs involving GW signals modeled after the circular TaylorEt approximant and the circular 
templates based on the TaylorT4 approximant, both requiring 2.5PN-accurate reactive dynamics. The TaylorEt GW signals are 
constructed using the 3PN-accurate expression for d<f>/dt and the 2.5PN-accurate expression for d^/dt, while the the TaylorT4 
approximant requires the 2.5PN-accurate expression for dx/dt. These PN-accurate expressions are extractable from Eqs. (6) 
and (7) in Ref. [33|. We conclude that for low mass binaries, the TaylorT4 approximant is not very effectual and in the case of 
high m binaries, the circular templates are fairly unfaithful. 



771l : 7712 


FF 


mt 


Vt 


5-5 


0.933 


9.900 


0.250 


3-9 


0.868 


9.888 


0.248 


3-12 


0.897 


11.26 


0.249 


5-10: 


0.903 


13.47 


0.250 


5-15: 


0.937 


16.20 


0.249 


10-10 


0.995 


20.51 


0.232 


15-15 


0.995 


30.72 


0.231 



TABLE V: The results originating from one of our FF experiments motivated by the Fig. 4 in Ref. [2(|. The inspiral 
GW signals are based on the circular TaylorEt approximant, given by Eq. (7) in Ref. [33| and therefore specified by the 
3PN-accurate expression for dcj)/dt and the 3.5PN-accurate expression for d£/dt. The circular templates are given by the 
2PN-accurate TaylorT4 approximant. Somewhat higher FFs are observed only for high mass binaries with m > 30 Mq. 



7711 /Af© : 7712 /Mq 


FF 


m t 


Vt 


5-5 


0.849 


9.744 


0.247 


1.4-10 


0.877 


6.891 


0.250 


3-9 


0.868 


9.871 


0.248 


3-12 


0.858 


14.058 


0.163 


10-10 


0.942 


18.940 


0.247 


15-15 


0.971 


27.701 


0.241 



TABLE VI: The initial LIGO FF values in the mass range relevant for the GRB 070201 progenitor following Ref. ||. The 
fiducial inspiral GW signals are from compact binaries having initial et = 10 -3 and computed using the phasing formalism of 
Sec. HI] The value of et can be treated to be representative of still lower residual initial eccentricities and the circular templates 
are based on the 2PN-accurate TaylorT4 approximant. Even though, the employed circular templates are neither effectual 
nor faithful for several mass combinations, we reiterate that further investigations involving the relevant LIGO data should be 
useful to make any physical implications of our ideal FF experiments. 



mi/M© 


7712 /Mq 


FF 


mt 


Vt 


3- 1 




0.837 


3.363 


0.250 


1-10 




0.909 


5.702 


0.250 


3-10 




0.876 


10.335 


0.249 


1-20 




0.937 


7.634 


0.250 


3-20 




0.945 


14.029 


0.250 


1-40 




0.966 


11.680 


0.195 


3-40 




0.972 


20.675 


0.208 



